Nonlinear elasticity of composite networks of stiff biopolymers with flexible linkers 
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Motivated by recent experiments showing nonlinear elasticity of in vitro networks of the biopoly- 
mer actin cross-linked with filamin, we present an effective medium theory of flexibly cross-linked 
stiff polymer networks. We model such networks by randomly oriented elastic rods connected by 
flexible connectors to a surrounding elastic continuum, which self-consistently represents the behav- 
ior of the rest of the network. This model yields a cross-over from a linear elastic regime to a highly 
nonlinear elastic regime that stiffens in a way quantitatively consistent with experiment. 
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The mechanical response of living cells depends largely 
on their cytoskeleton, a network of stiff protein polymers 
such as F-actin, along with various associated proteins for 
cross-linking and force generation. In addition to their 
importance for cell mechanics, cytoskeletal networks have 
also demonstrated novel elastic properties, especially in 
numerous in vitro studies 0, 0, 0, 0, 0, • The cellular 
cytoskeleton, however, is an inherently composite struc- 
ture, consisting of elements with highly varied mechan- 
ical properties, and there have been few theoretical or 
experimental studies of this aspect 1,0,13 El, El. Re- 
cent experiments on F-actin with the physiological cross- 
linker filamin have demonstrated several striking features 
while their linear modulus is significantly lower than for 
rigidly cross-linked actin systems, they can nonetheless 
withstand remarkably large stresses and can stiffen by a 
factor of 1000 with applied shear @, E3, El- This be " 
havior appears to result from the highly flexible nature 
of filamin, although the basic physics of such a network, 
in which the elasticity is dominated by cross-linkers, is 
not understood. Apart from their physiological impor- 
tance, such networks suggest new principles that may be 
extended to new synthetic materials with designed cross- 
links 0. 

Here, we develop a theoretical model for composite 
networks of rigid filaments connected by flexible cross- 
linkers, in which the macroscopic network elasticity is 
governed by the cross-links. We examine this model 
in a limit in which the basic elastic element is a single 
rigid rod, directly linked by numerous compliant cross- 
linkers to a surrounding linear elastic medium. We show 
that such a network stiffens in a manner determined by 
the mechanics of individual cross-links, which we model 
both as linear springs with finite extension, and also as 
wormlikc chains. We analyze our model in both a fully 
3D network, as well as a simplified ID representation, 
which already captures the essential physics of the non- 
linear behavior. The finite extension Iq of the cross-links 
along with the length of the filaments/rods L implies 
that there exists a characteristic strain 7 C ~ Eq/L for 




FIG. 1: (Color online) Schematic figure of an isotropic stiff 
polymer network with highly compliant cross- linkers. The in- 
set illustrates the proposed non-uniform deformation of the 
cross-linkers on a single filament in a sheared background 
medium. 



the onset of the nonlinear response of the network. In- 
deed prior in vitro experiments, in which the length of 
the cross-linkers was varied Q , have reported this linear 
dependence on Iq. We extend this model in a fully self- 
consistent manner, replacing the embedding medium by 
an effective medium whose elastic properties are deter- 
mined by those of the constituent rods and linkers. This 
self-consistent model can quantitatively account for the 
nonlinear response found in prior experiments on actin- 
filamin networks [1, [l(| • 

In a flexibly cross-linked stiff polymer network, ran- 
domly oriented stiff filaments or rods are intercon- 
nected by relatively short but highly flexible cross-linkers 
(Fig. [1} . The compliance of this network is dominated by 
the flexible cross-linkers, while the much stiffer filaments 
act mainly as a scaffold for the cross-linkers, ensuring 
rigidity of the network as a whole. Recent experiments 
have demonstrated that flexible biological cross-linkers 
such as filamin can be described as a semiflexible poly- 
mer using the wormlike chain (WLC) model 17 , 18| . The 
cross-linkers are characterized by their contour length £q 
and persistence length l v [l5| . A realistic force-extension 
curve of a typical biological cross-linker is shown as a 
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black solid line in the inset of Fig. [2j It is instructive 
to simplify this curve by assuming linear response with 
a spring constant fc c ; and a finite extensibility £ . This 
simplification retains the essential features, and is shown 
in Fig. [5] as a blue dashed line. We refer to this as simple 
cross-link behavior. 

To determine the elasticity of the network we use an 
effective medium approach, and divide the network into 
two mechanically connected sub-systems. The first con- 
sists of a stiff filament of length L decorated by n flexible 
cross-linkers, which we refer to as a hairy rod (HR). The 
other is the network connected to it, which we treat as 
an elastic continuum. Although the medium is assumed 
to deform affinely, we allow the local strain of the cross- 
linkers to depend on their position on the HR. By averag- 
ing over all orientations we may express the macroscopic 
stress in terms of the tension in a single HR connected 
to a medium, which is subject to a ID strain e along its 
backbone. The tension r in the center of this HR is the 
sum of the forces exerted by all cross-linkers on one half 
of the rod. To calculate these forces we treat the cross- 
linker as a spring connected in series with the medium, 
which we describe with a spring constant Kem- We are 
primarily interested in densely cross-linked networks for 
which Kem ~S> k c i- The extension of the cross-linker- 
medium system is given by ex at a distance x from the 
center of the rod. If the cross-linkers are homogenously 
distributed over the rod with a high density n/L we can 
write the sum over forces as an integral, 
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For strains e < j% only the first integral is present and 
the integration extends to L/2. In this case, the tension 
depends linearly on the strain. Using Eq. we compute 
the ID modulus G±d — T~/e, which is shown as a dashed 
blue line in Fig. O For small strains the system is linearly 
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strain j% a cross-over occurs to a second linear regime 

in which G\n asymptotically approaches ^uKemL. 

The nonlinear response of a cross-linker is more real- 
istically modeled with the WLC model (Fig.©. We 
calculate the tension in a rod with WLC cross-linkers 
analogously to Eq. ([I]). The ID modulus G\d is shown 
as a purple dash-dotted line in Fig. [2l Though quite 
similar to the simple cross-linker model, the more real- 
istic force-extension curve has introduced a considerable 
smoothing of the cross-over resulting in a gradual onset 
of nonlinear behavior of the HR with WLC cross- linkers. 
Nevertheless, the characteristic strain e c for the nonlin- 
ear behavior is proportional to £q/L independent of the 
exact nonlinear response of the linkers. 




FIG. 2: (Color online) The shear modulus G normalized by 
the linear modulus Go as a function of strain e (ID) or 7 (3D) 
for simple cross-linkers and WLC cross-linkers in the ID and 
the 3D version of the linear medium hairy rod model. In this 
plot we have chosen Kem = 100fc c ; as an example. The upper 
left inset shows the force extension curve of a simple cross- 
linker (blue dashed curve) and of a WLC cross-linker (black 
solid curve) . The lower right inset shows the normalized shear 
modulus as a function of strain 7 for various ratios of L/£q 
calculated with the self-consistent model. 



Using the ID model presented above we can compute 
the macroscopic stress of a network. A 3D isotropic net- 
work with a polymer length density p is modeled by an 
effective medium consisting of randomly oriented HR's. 
We can compute the macroscopic stress cr and shear mod- 
ulus G = 17/7 by averaging over all orientations 0, [||. 
The shear modulus is shown in Fig.[2]for the simple cross- 
linkers and for the WLC cross-linkers. The 3D curves 
are largely similar to the ID results, save for a factor two 
shift which may be understood by noting that the rods 
at a 45° angle to the stress plane, which bear most of the 
stress, experience an extensional strain e of 7/2. 

At large strains, when many of the cross-linkers are 
extended well into their nonlinear regimes, it is no longer 
realistic to assume a linear background medium. To ad- 
dress this, we shall now require the elasticity of the back- 
ground medium to self- consistently represent the nonlin- 
ear elasticity of its constituent HR's. 

Under strain the cross-links deform the surrounding 
elastic medium. The resulting longitudinal displacement 
S£ of the medium leads to a restoring force per unit length 
along the rod given approximately by the shear stiffness 
^ x 6£ The shear modulus of the medium depends 

on the density p of rods and the longitudinal stiffness 4^ , 
where e is the ID extensional strain of the medium along 
the rod: 



da ^ dr 

g?7 de 



(2) 



Here, A is a dimensionless geometric factor that depends 
on the architecture of the network. For an isotropic net- 
work in 3D, this is 1/15. (We note that this is a small- 
strain approximation and that A will be different for 
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anisotropic networks.) Thus, the effective stiffness Kem 
per cross-link is given by 



Kem = Ap- 



Ldr 

'de' 



(3) 



When subject to a shear strain 7, the resulting stress a 
within a network of rods can be expressed in terms of the 
tension r in each rod, which depends on its orientation 
relative to the shear plane. It is given by 



r(e) = 



L 



dx : 



de 1 



(4) 



where k c i (81) is the derivative of the force-extension curve 
of the cross-linker. Equivalently, we may write for r(e) 



dr 
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We first investigate the properties of this model using 
the simple force-extension curve (see inset Fig. [2]). For 
a densely cross-linked network we find a linear regime 
below 7 C = jfj2' F° r l ar S er strains the system enters a 
highly nonlinear regime for which 



dr 
de 
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This is in marked contrast with the linear medium model 
in which there is only a cross-over between two distinct 
linear regimes. 

A real network with compliant cross-linkers is more 
realistically modeled by solving Eq. ([5]) (numerically) us- 
ing the WLC force-extension curve for the cross- linkers. 
The shear modulus in this case, computed exactly as be- 
fore, is graphed in the lower right inset of Fig. [51 At low 
strains G ~ nk c ipL and there is a gradual onset of non- 
linear response originating from the nonlinear entropic 
elasticity of the cross-linkers(see upper left inset Fig. 
At a strain ~ Iq]L the cross-linkers at the edges of the 
rods become effectively rigid, which marks the onset of 
the nonlinear network behavior. 

In view of the nonlinearity of this system, it is more 
appropriate to use a differential modulus K = ^ rather 
than G. The differential modulus is plotted in Fig. [3] 
Up to a critical stress o~ Cl the elasticity is dominated by 
WLC cross-linkers placed on a rigid rod connected to a 
much stiffer medium. At larger stresses, the cross-linkers 
at the edges of the HR reach full extension and, conse- 
quently, couple strongly to the surrounding network. In 
this limit, the slope in a log (if) vs log(cr) plot approaches 
«1- l/(jApL 2 — 1), as it does for simple cross-linkers. 
This exponent is a consequence of the composite nature 
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FIG. 3: (Color online) Differential modulus K = dr/dy nor- 
malized by the linear modulus Go as a function of stress a 
normalized by a c for the self-consistent (Self.-Con.) model. 
We also plot K/Go for the linear medium (Lin. Med.) model 
and a model for rigidly cross-linked semiflexible polymer net- 
works (Rig. C.-L.). The black line indicates a slope of 1. 
The inset shows the reduced tension profile (j> along the rod, 
normalized by the mid-point tension r in Eqs. (fl} and Jjj). 



of the network and its nonlinear constituents, although 
it is independent of the exact form of the nonlinear re- 
sponse of the cross-linkers. For a dense flexibly cross- 
linked network pi? ^> 1 and, therefore, we expect a slope 
of 1. This is consistent with recent experimental data 
on actin networks cross-linked by the highly compliant 
cross-linker filamin in which a slope of 1 was found Q in 
contrast to a slope of 3/2 found for rigidly cross- linked 
networks 0, Ll3| • Interestingly, in vivo experiments show 
that cells also exhibit powerlaw stiffening with an expo- 



nent of 1 

We compare our results to the linear medium model 
and a model based on the nonlinear response of the semi- 
flexible actin segments between cross-links that has been 
used successfully to describe rigidly cross-linked actin 
networks @, E3| in Fig. H Although the three curves co- 
incide for small stresses, at intermediate stresses a > a c 
the linear medium model curve rolls over to a linear 
regime. Clearly, our self-consistent model and the model 
for rigidly linked networks begin to differ in the nonlinear 
regime. 

So far, we have considered only the mid-point tension. 
In networks of elastic filaments of finite length, however, 
the tension along a single filament is not uniform, but 
decreases towards its ends 



2J,|23J, 24]. The inset of Fig. 



[3] shows the ratio </> of the tension at point x along the 
rod to the maximum tension. This maximum occurs at 
the mid-point x = 0, and is given by r in Eqs. (|4]) and 
([5]). The tension at a point x can be obtained by re- 
placing the lower limit of of the x' integral in Eq. ((4]) by 
x. The tension profile is parabolic below j c and quickly 
converges to a more flattened out profile in the nonlinear 
regime. We can use the tension profile to relate the max- 
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imum tension in a single HR to the macroscopic stress 
a 12 1( . For typical experimental conditions in an actin- 



filamin gel[l3| we estimate a maximum force on a single 
cross-link to be at most 5 pN for isotropic rods and of 
order 1 pN or less for oriented rods. 

A feature shared by the linear medium model and the 
self-consistent model is the characteristic strain j c ~ 
4£q/L for the onset of nonlinear response. The propor- 
tionality with £ is consistent with the results of Wag- 
ner et al. where cross-linker length was varied, although 
they observed larger values of 7 C than expected either 
from our model and based on Refs.jH, [l(J. Recent ex- 
periments on actin-filamin networks also show a depen- 
dence of the critical strain that is approximately inverse 
in actin filament length L 14| in agreement with our re- 
sults. This sensitivity of network response to filament 
length, both in experiments and in our model, appears 
to be one of the hallmarks of actin-filamin networks. On 
the one hand, this may explain the apparent difference 
between the critical strains reported in Refs. [H, 0, fioj ] . 
On the other hand, it also suggests that it may be more 
important to directly measure the filament length dis- 
tribution in such experiments than in other similar in 
vitro studies. In Wagner et al., for instance, the filament 
length was not measured, but was inferred from prior re- 
ports of the length dependence on the capping protein 
gelsolin 0]. 

In previous work, DiDonna and Levine have assumed 
a sawtooth force-extension curve for the cross-linkers to 
mimic domain unfolding. They report a fragile state with 
shear softening when an appreciable number of cross- 
linkers are at the threshold of domain unfolding llj . Our 
model is based on the stiffening of the cross-linkers, which 
occurs at forces far below those required for domain un- 
folding [13, EH • This leads to strain stiffening at a point 
where only a fraction of cross-linkers are at their thresh- 
old for nonlinear response. Thus in both our model and 
that of Ref. the network responds strongly to small 
strain changes, though in an opposite manner: stiffen- 
ing in the present case vs. softening in Ref. llj. In 
related work, Dalheimer et al. show that isotropic net- 
works linked by large compliant cross-linkers exhibit a 
shear induced ordering transition to a nematic phase 
Our model accounts for the architecture of the network 
through an averaging procedure in a scalar quantity A. 
We are presently investigating the effect of an ordering 
transition on the nonlinear response of the network. 

We have introduced a model for flexibly cross-linked 
stiff polymer networks based on cross-linker elasticity. 
Our model yields an exponent of 1 in the asymptotic 
powerlaw behavior of a K vs a curve in agreement with 



experiments on in vitro filamin-actin networks @, [Io| . 
The exact form of the nonlinear response predicted by 
our model can be tested by further experiments 
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